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Membrane targeting domains play crucial roles in the association of signalling molecules to the 
plasma membrane. For most peripheral proteins, the protein-to-membrane interaction is transient. 
After proteins dissociate from the membrane they have been observed to rebind following brief 
excursions in the bulk solution. Such membrane hops can have broad implications for the effi¬ 
ciency of reactions on membranes. We study the diffusion of membrane-targeting C2 domains using 
single-molecule tracking in supported lipid bilayers. The ensemble-averaged mean square displace¬ 
ment (MSD) exhibits superdiffusive behaviour. However, traditional time-averaged MSD analysis 
of individual trajectories remains linear and it does not reveal superdiffusion. Our observations are 
explained in terms of bulk excursions that introduce jumps with a heavy-tail distribution. These 
hopping events allow proteins to explore large areas in a short time. The experimental results are 
shown to be consistent with analytical models of bulk-mediated diffusion and numerical simulations. 


A myriad of signalling proteins are recruited to specific 
cell membranes via phospholipid-binding domains [HE]. 
These molecules dock to the surface of specific lipid mem¬ 
branes and undergo two-dimensional diffusion in search 
of a target. Once the target is located, many proteins 
either activate or suppress a downstream signalling path¬ 
way for various physiological and pathological processes. 
Examples of membrane-targeting domains include pleck- 
strin homology (PH) [3 and C2 [4], which have been 
identified in hundreds of human signalling molecules as 
well as in eukaryotic species as diverse as fungi and flies 
[Q. PH domains bind specifically to phosphoinositides 
while C2 domains bind a variety of membranes, and a 
subset of C2 domains only bind membranes in the pres¬ 
ence of calcium and play key roles in signalling pathways. 
The association to lipid membranes often takes place in 
response to different extracellular and intracellular stim¬ 
uli, but typically the residence on the membrane surface 
is only temporary. The transient nature of peripheral 
protein-membrane interactions enables a tight temporal 
regulation of signal transduction. Further, membrane 
dissociation has also broad implications on the search for 
the target substrate, but this process is less understood. 

Recently, Knight and Falke observed the dissociation 
of PH domains from supported bilayers followed by rapid 
rebinding to the surface after a short excursion in the 
bulk solution [B]. They proposed that the hopping pro¬ 
cess may be important in the search for target molecules 
in eukaryotic cells. Subsequently, Yasui et al. found that 
PTEN (phosphatase and tensin homologue) molecules 
hop along the plasma membrane of living cells due to 
dissociation followed by rebinding [7j. PTEN is an im¬ 
portant protein that suppresses development of cancer. 
It prevents cells from growing and dividing too rapidly 
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by dephosphorylating phosphoinositide substrates on the 
plasma membrane. PTEN-membrane affinity is regu¬ 
lated by a C2 domain and it is enhanced by electrostatic 
interactions. The observed hopping of the C2 domain 
on the plasma membrane is thus expected to alter the 
dynamics of the search for a phospholipid substrate. 

A straightforward consequence of membrane hopping 
is that a molecule remains in its immediate vicinity for 
a short time and then jumps to a location that is fur¬ 
ther away than expected from two-dimensional diffu¬ 
sion. Therefore, the search process is allowed to explore 
larger areas and the molecule can bypass diffusion bar¬ 
riers that may be present in the membrane. However, 
hopping comes at the cost of the search being less ex¬ 
haustive. We may ask the questions how the dynamics 
of membrane-targeting domains is affected by such long 
jumps and how this motion deviates from a simpler two- 
dimensional diffusion. Such potential complex behaviour 
can yield anomalous diffusion of membrane-targeting do¬ 
mains, which would alter the outcome of search processes 
and the sequential molecular reactions. 

Anomalous diffusion is widespread in the motion of 
molecules in biological systems [MU- In general, a par¬ 
ticle exhibits anomalous diffusion when the mean square 
displacement (MSD) scales as a power law with an expo¬ 
nent a/1 

(x 2 (t)) = K a t a , (1) 

where K a is the generalized diffusion coefficient with 
units cm 2 /s“. When a < 1 the process is subdiffusive 
and when a > 1 it is superdiffusive. Subdiffusion in 
the cytoplasm jl2Hl4] . the nucleus US], and the plasma 
membrane jlBHTS] of live cells is caused by crowding 
mm and complex interactions with the cytoskeleton 
and macromolecular complexes, among others. Similarly, 
subdiffusion can take place in model membranes due to 
crowding and packing effects mm- The appearance 
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FIG. 1. Sketch of the diffusion process. A molecule alternates 
between phases of two-dimensional and three-dimensional dif¬ 
fusion. Diffusion in the three-dimensional bulk is much faster 
than diffusion on the lipid bilayer, and thus only the effective 
two-dimensional process is observed without loss of trajectory 
connectivity. The excursions into the bulk are seen as long 
jumps in the two-dimensional trajectories. 


of superdiffusion processes in biomolecular systems is far 
less common. The archetypal mode of superdiffusive mo¬ 
tion is due to active cytoplasmic flows and transport me¬ 
diated by molecular motors, requiring ATP energy con¬ 
sumption [25H25| ■ 

From a theoretical point, there are three major mech¬ 
anisms that can introduce superdiffusion |26) . It can 
be caused by correlations in the random walk, such as 
those in fractional Brownian motion with a Hurst in¬ 
dex H > 1/2, by persistent directional motions (Levy 
walks), and by long jumps (Levy flights). Active biolog¬ 
ical transport can be modelled as Levy walks [23] • Bulk- 
mediated diffusion processes, which can be described as 
Levy flights, have been observed for transient adsorption 
on a solid surface where molecules display intermittent 
behaviour, alternating between periods of immobilization 
at the solid-liquid interface and periods of diffusion in the 
bulk fluid [27j [2|]. 

In this article we report the experimental observation 
of superdiffusive transport of membrane-targeting C2 do¬ 
mains on supported lipid bilayers. Measurements of the 
diffusion of membrane-targeting domains are performed 
by single-particle tracking and are compared to both an¬ 
alytical theory and numerical simulations. In stark con¬ 
trast to active cytoplasmic transport, superdiffusion in 
model membranes does not require energy. Our data 
strongly suggests that superdiffusion is caused by bulk- 
mediated diffusion, namely molecules dissociate from the 
membrane and perform three-dimensional random walks 
until they reach the membrane again and readsorb at 
a new location, as sketched in Figure [l] Interestingly, 
the motion of membrane-targeting domains shows weak 
ergodicity breaking, a phenomenon that has recently at¬ 
tracted considerable attention in cellular environments 
and other complex systems 0(11 (MSB- The ergodic 
hypothesis, which is fundamental to statistical mechan¬ 
ics, states that ensemble averages and long-time averages 
of individual trajectories are equivalent. The violation of 
ergodicity has pronounced implications for the dynamics 


of individual molecules, which can be very different from 
the ensemble statistics [8]. In the traditional way of ob¬ 
taining the MSD, the square displacements are averaged 
over a large ensemble of molecules at a time t since the 
beginning of the measurement, i.e. an ensemble average. 
Alternatively, it is possible to perform the average over 
all the displacements in a lag time A of a single trajec¬ 
tory, i.e. a temporal average. For ergodic systems, both 
averages converge to the same value. However, weak er¬ 
godicity breaking can take place as a consequence of ki¬ 
netics with power-law statistics in the plasma membrane 
[521 [55] and in the cytoplasm of live cells [HIM] as well 
as in inorganic complex systems such as quantum dots 
[35l ;36j and models of glassy dynamics \29\. 


RESULTS 

Diffusion of membrane targeting proteins on 
supported lipid bilayers 

We tracked the motion of the membrane-targeting 
C2A domain from synaptotagmin 7 [57], labelled with 
Atto-565, in a supported lipid bilayer composed of phos¬ 
phatidylcholine (PC) and phosphatidylserine (PS) at a 
3:1 ratio. The lipid bilayer was self-assembled on a 
clean coverslip 0. Imaging was done in a home-built 
total internal reflection (TIRF) microscope under contin¬ 
uous illumination at 20 frames/s. Surface densities were 
kept low enough to enable accurate tracing of trajecto¬ 
ries and to allow assignment of connections even after 
micrometer-long jumps. 

Figure [2]i shows an example of trajectories obtained in 
a 10-s window, overlaid on the last frame. Often, long 
jumps are observed in the particle trajectories as seen in 
the examples in Figs. and c. These jumps suggest 
the C2A molecules detach from the surface and readsorb 
after brief excursions into the liquid bulk. The motion 
in the bulk is much faster than diffusion on the viscous 
membrane and jumps are thus expected to occur instan¬ 
taneously for all practical purposes. For the C2A domain, 
the diffusion coefficient in the lipid bilayer D s is of the 
order of 2 /im 2 /s, but in liquid the diffusion coefficient 
Db is estimated to be 100 times higher 138]. As a con¬ 
sequence, when a molecule performs a jump through the 
bulk it can sometimes be observed at reduced intensity in 
both the old and new locations within the same imaging 
frame, as seen in Figure [2|:>. 

In order to study the effect of the dissociation con¬ 
stant, we also employed a C2A construct fused to a non- 
membrane interacting glutathione S-transferase (GST), 
which has a strong tendency to dimerize (Figure |3^i). 
The GST-C2A dimer forms two independent interactions 
with the membrane and will consequently have a slower 
dissociation rate than C2A monomer, providing a good 
comparison for validating our superdiffusion predictions. 
Additionally, GST-C2A dimer has a higher viscous drag 
coefficient and, in turn, its diffusion coefficient on the 
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FIG. 2. Single particle tracking of membrane-targeting do¬ 
mains. (a) C2A-Atto565 individual trajectories collected dur¬ 
ing a 10-s time window. Three mobile trajectories are ob¬ 
served in the image together with one immobile particle that 
is tracked but is not included in the analysis. The last frame 
is overlaid on the trajectories. The original data is shown in 
Supplementary Video SI. Scale bar 2 pm. (b) Region of in¬ 
terest (ROI) around the location of a micrometer jump that 
occur in the lowermost trajectory, marked with the letter b. 
Three frames are shown corresponding to before, during, and 
after the jump. Scale bar 0.5 pm. (c) ROI around the location 
of the jump marked with the letter c. Scale bar 1 pm. 


membrane surface is reduced to nearly half [39]. 

We collected 14,000 C2A and 3,600 GST-C2A mobile 
trajectories. Immobile fluorophores that did not exhibit 
any apparent diffusive motion were excluded from the 
analysis. The ensemble-averaged MSD ( r 2 (t )) of C2A 
monomers and dimer-forming GST-C2A are shown in 
Figure A deviation from a linear MSD is evident 
in the figure, showing superdiffusive behaviour. Further, 
the onset of superdiffusion for GST-C2A occurs at a later 
stage. 

The time-averaged MSD S 2 ( A) is often used in the 
analysis of individual trajectories. Throughout this 
manuscript we will denote the ensemble average of an 
observable with brackets (•) and the time average with 
an overbar 7 . For a trajectory with N time points, 

N-n 

<5 2 ( A ) = [ r C? V + A ) - r 0V)] 2 , (2) 

i=i 

where r is the time interval between consecutive mea¬ 
surements and n = A/r. This approach is especially 
useful when a limited number of trajectories is available, 
as usually occurs in single-molecule studies. Figure 0 



b 


CM 

E 



CM 

E 

3 

1— 

a. 



FIG. 3. Anomalous diffusion analysis of membrane-targeting 
domain C2A (monomer) and dimer forming GST-C2A. (a) 
Sketch of the C2A monomer and the GST-C2A dimer em¬ 
ployed in this study, (b) Ensemble averaged MSD (r 2 (t)). (c) 
Time averaged MSD (<5 2 ) as a function of lag time A. The 
time average-MSD of individual trajectories varies greatly, so 
the MSDs of individual trajectories are also ensemble aver¬ 
aged. (d-e) Distribution of displacements for A = 100 ms. 
The total number of displacements are 207,000 and 56,000 for 
C2A and GST-C2A, respectively. The solid lines show fitting 
to equation 0 and to the individual components of the propa¬ 
gator, i.e. the Gaussian part [(1 —w)/27rcr 2 ]exp(—r 2 /2cr 2 ) and 
the Cauchy propagator part ury/27r(r 2 -|-y 2 ) 3 / 2 . The cutoff at 
2.6 pm appears because trajectories are not connected when 
jumps longer than this distance take place. This threshold is 
placed in order to avoid the risk of particle misconnections. 


shows the time-averaged MSD after it is additionally av¬ 
eraged over all the trajectories. GST-C2A exhibits the 
expected slower diffusion rate than C2A, based on the 
MSD slope. As mentioned above, for ergodic processes, 
the temporal and ensemble averages coincide in the long 
time limit, <5 2 (A) = (r 2 (A)). However, the ergodic hy¬ 
pothesis breaks down for C2A molecules. In contrast to 
the ensemble-averaged MSD, the time-averaged MSD is 
linear in lag time 

(TO) ~ A - (3) 

Thus, an observer analysing time-averages would reach 
the misleading conclusion that the diffusion behaviour is 
not anomalous. 

The distribution of displacements P(r) at A = 100 ms 
is shown in Figs. [3ji and e for C2A and GST-C2A, respec- 
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tively. The distribution exhibits two different character¬ 
istic regimes: a central part up to a distance r ss 1.5 /im 
and a long tail. This behaviour can be understood 
from the scaling properties of bulk-mediated diffusion 
as discussed by Bychuk and O’Shaughnessy m- Once 
a molecule dissociates from the surface, it performs a 
three-dimensional random walk until it returns. In the 
asymptotic limit, the first return time distribution scales 
as ~ T” 1,5 . For any given return time, the surface 
distance between the dissociation and return points has a 
Gaussian distribution P(rj\r) ~ exp(—rJ/4D & r). There¬ 
fore, the distribution of jump lengths is P( 17 ) ~ rj 3 , as 
observed in Figs. P and e for long distances. 

The theoretical probability density function of jump 
lengths can be found using the image method [JT]. The 
distance of first return to the surface are governed by 
-P(r) = 7 o/ 2 tt (r 2 + 7 q) 3/2 , that is a two-dimensional 
Cauchy distribution. At short times, the probability that 
the particle performs more than a single jump is small. 
If we neglect the distance covered by surface diffusion 
within time intervals at which the particle undergoes a 
bulk excursion, the motion at each short interval is ei¬ 
ther by surface diffusion or via a jump. We can then 
approximate the distribution of displacements at short 
times by 


P( r) 


7o 

27r(r 2 + 7g) 3 / 2 


(! ~v) 

27T<7 2 


exp(—r 2 / 2 tr 2 ), 


(4) 
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FIG. 4. Numerical simulations of Levy flights. 500 realiza¬ 
tions were performed, in which a particle alternates between 
2D random walks and bulk-mediated jumps, (a) Probability 
density of the tracer displacements. The density is well de¬ 
scribed by a theoretical model that includes a Gaussian cen¬ 
tral part and a Cauchy propagator of the form 7 o/(r 2 + 7 o) 3 ^ 2 - 
(b) Ensemble-averaged MSD (r 2 (t)} as a function of time. 
The ensemble-averaged MSD is computed from the distance 
covered by the tracer in a time t from the start of the realiza¬ 
tion. (c) The time-averaged MSD <5 2 (A) is averaged over all 
realizations and plot against lag time A. 


where u> is the probability that the particle hops within 
the given time and surface diffusion yields er 2 = 2 D s t. A 
least-square fitting of the distribution of displacements 
(Figs.p and e) to this propagator yields D s = 1.7 /im 2 /s 
for C2A monomers and D s = 1.0 pm 2 /s for GST-C2A. 
The parameter 7 is found to be 0.24 /rrn and 0.12 /im for 
C2A and GST-C2A, respectively. 

The distribution of displacements for longer times in¬ 
volves both a random number of jumps, each having a 
Cauchy distribution, and the Brownian motion on the 
surface. Chechkin et al. derived the full solution for the 
propagator of bulk-mediated diffusion [42] . For the case 
when D s = 0 and neglecting long distance corrections, 
the distribution of displacements is given by the Cauchy 
propagator, in agreement with scaling arguments [401 . 

P ^ ~ 27rlr 2 ^(7tfj 3 ^' ^ 

When the particles also diffuse on the surface, i.e. D s ^ 0 
the probability density of the displacements is given by 
the convolution of equation © with a normal distribu¬ 
tion. Even though the full solution for long times is com¬ 
plicated, the tail of this distribution for large distances 
still scales as P( r) ~ r~ 3 . Due to this asymptotic be¬ 
haviour, the exact distribution has similar properties to 
the Cauchy distribution. 


Numerical simulations: diffusion in the presence of 
bulk excursions 

In order to verify the model of surface diffusion in the 
presence of bulk excursions we analyse numerical simula¬ 
tions of the process diagrammed in Figure [I] Molecules 
perform a two-dimensional random walk, but at random 
times they jump due to a hypothetical bulk excursion. 
The surface residence times are assumed to be indepen¬ 
dent and identically distributed exponential random vari¬ 
ables and the jumps are modelled according to the first 
return time to the surface given simple diffusion in a 
three-dimensional medium. These simulations are anal¬ 
ysed in the same way as with experimental observations 
of the motion of membrane-targeting C2 domains on sup¬ 
ported membranes. 500 realizations were simulated off- 
lattice with a surface diffusion coefficient D s = 0.5 and 
a dissociation coefficient k = 0.1. The chosen param¬ 
eters do not intend to capture the real protein prop¬ 
erties, but to simply test theoretical predictions with¬ 
out the effects of experimental noise. The displacements 
for two-dimensional diffusion are drawn from a Gaus¬ 
sian distribution with variance cr 2 = 1 and the return 
times from bulk excursions are drawn from a distribution 
ip(tb) = zo( 47 rDbf 3 ) - V 2 exp(- 2 o/ 4 D b f b ) [41] . Then the 
jump distances are drawn from a Gaussian distribution 
with variance er 2 = 2 D b t b . 
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The distribution of displacements P(r) for the numeri¬ 
cal simulations is shown in Figure [4^,. As expected, there 
are two regimes: a central Gaussian part due to the two- 
dimensional diffusion on the membrane between bulk ex¬ 
cursions, and a heavy tail that arises from the long dis¬ 
tance behaviour of bulk excursions. The distributions 
for short times can again be modelled with a propagator 
that includes contributions from Gaussian surface diffu¬ 
sion and a Cauchy distribution due to bulk excursions. 
By fitting to equation Q, it is found D s = 0.50 ± 0.05 
(the value employed in the simulations is D s = 0.5) and 
70 = 0.75. 


MSD analysis 

The dynamics of a particle with a Cauchy propagator 
are particularly interesting because the theoretical vari¬ 
ance of the displacements diverges, 

nOO 

(r 2 (f)) = / (2nr)r 2 P(r)dr = oo. ( 6 ) 

Jo 

In practice, a diverging second moment implies that there 
is a non-negligible probability for the occurrence of ex¬ 
tremely long jumps and this phenomenon has direct im¬ 
plications in the measured MSD. Figure shows the 
ensemble-averaged MSD computed from the numerical 
simulations. The MSD increases in a superlinear fash¬ 
ion, i.e. by employing equation ([l]), we have a > 1, 
which implies the process is superdiffusive. 

Let us now analyse the unexpected MSD behaviour, 
starting from the time-averaged MSD of individual tra¬ 
jectories. We can show that the time-averaged MSD is 
linear in lag time for any random walk with independent 
increments Uj = r, +1 — 17 , such that (u* • u ; ) = 0 when 
i ^ j. From the definition of the time averaged MSD 
(equation ([ 2 ])) (35] , 


d 2 (A) = 


where we have used the approximation that t > A, we 
omitted the term u * ' u j because it is zero on 

average, and again we have used the parameter n = A/r. 
Therefore we see that for symmetric random walks with 
independent increments, the time-averaged MSD is linear 
as observed in Figs. and 15 - 

Although the time-averaged MSD for individual tra¬ 
jectories is linear, the ensemble averaged MSD ( r 2 (t )) is 
not. We can understand the superdiffusive behaviour by 
assuming we can define the motion in terms of two in¬ 
dependent processes r (f) = b (t) + y (t), where b(f) is a 


N-n 

1 

N-n 

N-n 


A 
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N-n 

/ i+n \ 2 


E 

E u 0 

(7) 

2—0 

\k=i / 


N—n i 

+n 


E E u * 

( 8 ) 

2 — 0 / 

C=2 


u i, 


(9) 


i—1 


two-dimensional Brownian motion and y (t) is a Levy pro¬ 
cess with a probability density defined by equation ([5|. 
Then the MSD is (r 2 ) = (b 2 ) + (y 2 ). The first term is 
linear in time but the second term has a superdiffusive 
nature [213 45} j44] . 


DISCUSSION 

The propagator for surface diffusion in the presence 
of bulk-mediated jumps (equation Q) depends on the 
surface diffusion coefficient D s and the parameter 7 that 
reflects the transition between the surface and the bulk 
phase. Namely, 7 ~ a/r ^ es , where Td es is the mean 
desorption time and a is a dimensional factor. Bulk- 
mediated diffusion thus predicts 7dimer < Tmonomer, in 
agreement with the values we find for C2A and GST- 
C2A. 

The surface motion of these membrane-targeting do¬ 
mains is well described by Levy flights, a random walk 
where the step displacements have a heavy-tailed dis¬ 
tribution. The heavy tail arises from the dissociation 
of molecules from the membrane, which then perform 
a three-dimensional random walk until they reach the 
surface again at another location. The process involves 
the first return to a surface and it converges to a power 
law according to the Sparre-Andersen theorem j41j . This 
type of Levy flight dynamics is fundamentally different 
from Levy walks induced by molecular motors in the cy¬ 
toplasm because periods of active motion require an en¬ 
ergy input, typically in the form of ATP hydrolysis, while 
bulk excursions occur spontaneously. 

One of the most interesting effects of the observed 
bulk-mediated diffusion statistics is that the ensemble- 
averaged MSD exhibit superdiffusive behaviour, whereas 
the temporal averages suggest normal diffusion. This 
nonergodic behaviour is similar to that of continuous 
time random walks (CTRW) where the sojourn time dis¬ 
tribution between steps has a probability distribution 
that is heavy-tailed. Also in the CTRW, (<5 2 (A)) ~ A 
and ( r 2 (t )) ~ t a , albeit the CTRW is subdiffusive with 
a < 1. The difference in the behaviour of temporal and 
ensemble averages is the key signature of weak ergodicity 
breaking in the process |45j . 

To date, different groups have observed normal diffu¬ 
sion for membrane proteins in supported lipid bilayers, 
which appear to contradict our findings [4HH451 . There 
are several reasons for this apparent discrepancy. Single¬ 
particle tracking in lipid bilayers often focuses on time- 
averaged MSD, which does not show any non-linearity 
in lag time. Thus it would be reasonable to reach the 
conclusion that diffusion is not anomalous. Furthermore, 
anomalous diffusion in supported bilayers is known to de¬ 
velop as a result of packing and crowding. These mecha¬ 
nisms are modelled by a fractional Langevin equation, 
which is ergodic in nature, with anomalies that show 
up in the time averages. The distribution of displace¬ 
ments has also been previously reported as exhibiting 
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Gaussian behaviour. Here we report on the motion of 
surface-bound membrane domains that exhibit desorp¬ 
tion from the membrane within the experimental obser¬ 
vation time. The behaviour of transmembrane proteins 
or lipids is very different because the free energy barrier 
for release from the membrane is too high to be observed 
within the constrains of experimental observations @S1- 
148] . Previous works dealing with membrane-targeting 
domains such as C2 have generally been limited to short 
displacements in order to exclude the effect of long bulk- 
mediated jumps in diffusion measurements [551 H5] . 

What are the biological implications of surface su- 
perdiffusion for peripheral membrane proteins? Search 
processes are ubiquitous in cell biology and it is feasible 
to assume that evolution has optimized search param¬ 
eters. For signalling molecules delivered to the plasma 
membrane during a specific stimulus, the target molecule 
is often scarce in a sea of other lipids and proteins. Thus 
we can envision that if a molecule does not find its tar¬ 
get in a given time, it becomes more efficient to start 
searching at a different location. Is it appropriate then 
to assume Levy flights yield the optimal search for sparse 
targets when compared to Brownian motion? For one¬ 
dimensional intermittent processes that switch between 
Brownian motion and ballistic relocation phases, it has 
been shown that the search process is significantly more 
efficient when relocation times are power-law distributed, 
resulting in a Levy walk [50j . Notably, when Levy dy¬ 
namics are employed, the search is less sensitive to the 
target density EH- In general, the optimal strategy de¬ 
pends on the average target distance from the starting 
point m- However, blind searches inside a living cell are 
very different from a search in an unobstructed environ¬ 
ment. Several aspects provide additional complexities in 
the plasma membrane, in particular m- Experimental 
measurements show that the plasma membrane is com¬ 
partmentalized in a way that proteins and lipids have 
the tendency to remain transiently confined within small 
regions [52]. Further, membrane proteins typically ex¬ 
hibit subdiffusion with anti-persistent increments where 
molecules drift towards the locations that they visited 
in the past. While this subdiffusive behaviour provides 
the opportunity for a thorough and compact search, it is 
definitely not the optimum situation to find sparse tar¬ 
gets. A superdiffusive Levy flight provides a mechanism 
to overcome the effects of anti-persistent correlated sub- 
diffusive motion. Thus, we expect Levy flight dynamics 
to often outperform a Brownian search. 

The obstruction to the diffusion of membrane 
molecules has two different sources, both of them caus¬ 
ing anti-persistent correlations in the random walk. On 
one hand, obstacles can be introduced by immobile trans¬ 
membrane proteins which affect all lipids and membrane 
proteins. On the other hand, a more severe obstruction 
can be caused by cytoskeleton components that may not 
be in direct contact with the plasma membrane [53] . The 
effect of these barriers is not equal for all membrane pro¬ 
teins. Proteins that have large intracellular complexes 


are blocked much more efficiently than small molecules. 
In cases where a large signalling molecule adheres to the 
membrane via phospholipid-binding domains, bulk ex¬ 
cursions allow for the exploration of larger areas. Other¬ 
wise, the molecule would remain confined for long times 
within cytoskeleton-formed corrals, even when no sub¬ 
strate target is found within this region. 

In summary, we have observed the nonergodic, su¬ 
perdiffusive motion of membrane-targeting peptide do¬ 
mains in supported lipid bilayers. The motion is well- 
described by Levy flights with jumps that have a heavy- 
tail distribution. The long jumps are caused by excur¬ 
sions into the liquid bulk. After dissociating from the 
membrane, the molecules diffuse in three dimensions un¬ 
til they reach the membrane again and bind at a new 
location. Diffusion in the liquid bulk is much faster than 
diffusion in the membrane, therefore we do not consider 
the delay time between dissociation and readsorption. 
The surface distances covered by jumps have a Cauchy 
distribution, which is responsible for the heavy tail in the 
superdiffusive Levy flights. Model membranes provide an 
elegant system to study the effect of superdiffusive Levy 
flights because they are not subjected to the interactions 
with other cell components that would mask its experi¬ 
mental observation. However, hopping was already ob¬ 
served on the surface of live cells [7] and we foresee these 
processes have broad physiological relevance in the sur¬ 
face diffusion of signalling molecules. 

METHODS 
Imaging buffer 

Imaging and rinsing during the preparation steps was 
performed in an imaging buffer consisting of 50 mM 
HEPES, 75 mM NaCl, 1 mM MgCl 2 , 2 mM tris(2- 
carboxyethyljphosphine (TCEP), 200 /zM CaCl 2 . CaCl 2 
is necessary for C2 domain binding to the reconstituted 
membrane. 


Preparation of phospholipid vesicles 

Phospholipids were purchased from Avanti Polar 
Lipids (Alabaster, AL). Chloroform-suspended 18:1 (A9- 
Cis) PC (DOPC) and 18:1 PS (DOPS) were mixed at a 
ratio of 3:1. The phospholipid mixture was vacuum dried 
overnight and resuspended in imaging buffer to a final 
concentration of 3 mM followed by probe sonication to 
form sonicated unilamellar vesicles (SUVs). 

Preparation of coveslips and supported lipid bilayers 

Glass coverslips were cleaned by sonication in a de¬ 
tergent solution followed by soaking in 1M KOH. The 
coverslips were rinsed extensively in Milli-Q water and 
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blown dry with a stream of nitrogen gas. Then, the 
coverslips were treated with an oxygen plasma. Imme¬ 
diately after the plasma cleaning, a perfusion chamber 
(CoverWell, Grace Bio-Labs) was adhered to the cover- 
slip. In order to deposit the lipid bilayers a solution of 
SUVs (1.5-mM lipid) composed of phosphatidylcholine 
(PC) and phosphatidylserine (PS) at a 3:1 ratio in 1M 
NaCl and imaging buffer was introduced into the per¬ 
fusion chamber and incubated for one hour at 4°C. Re¬ 
frigeration minimizes lipid oxidation. The surface was 
then rinsed with imaging buffer multiple times prior to 
addition of protein sample. 

C2A and GST-C2A expression and purification 

An expression plasmid containing the gene for a GST- 
ybbR-Synaptotagmin 7 (Syt7) C2A domain fusion pro¬ 
tein was transformed into E. coli BL21-CodonPlus(DE3) 
competent cells. The ybbR segment provides a site for 
Sfp-catalysed fluorophore labelling [S3] . Cells were grown 
at 37°C to an OD 6 00 of 0.6 and then induced to express 
protein with 0.5 mM IPTG at room temperature for 6 
hours. The harvested cells were lysed at 18,000 lb/in 1 2 3 
in a microfhiidizer in a buffer containing 50 mM Tris pH 
7.5, 400 mM NaCl and centrifuged at 17,000 rpm in a 
Sorval SS-34 rotor. The clarified lysate was loaded onto 
a 5-rnl GSTrap FF column (GE Healthcare LifeSciences, 
Pittsburgh, PA) followed by gradient elution with 50 mM 
Tris, pH 8.0, 100 mM NaCl, and 10 mM glutathione. 
Fractions containing protein were pooled and diluted to 
reduce the salt to less than 0.1 M prior to loading onto 
a HiTrap Q HP column (GE Healthcare LifeSciences, 
Pittsburgh, PA) and eluting with a linear gradient to 
1 M NaCl in 25 mM Tris, pH 8.5, 20%(vol/vol) glycerol, 
and 0.02%(wt/vol) NaN 3 . A portion of the construct 
was subjected to thrombin cleavage and then separated 
using a Superdex 200 gel filtration column (GE Health¬ 
care LifeSciences, Pittsburgh, PA) equilibrated in 50 mM 
Tris, pH 7.5 and lOOmM NaCl to yield a ybbr-Syt7 C2A 
construct. 


Protein labeling 

10 mM CoASH (New England Biolabs, Ipswich, MA) 
in 400 mM Tris, pH 7.5 was mixed with 10 mM 


ATTO-565 maleimide (ATTO-TEC, Siegen, Germany) 
in dimethylformamide and incubated at 30° C overnight 
to form ATTO-565 CoA, then quenched with 5 mM 
DTT, 10 mM Tris pH 7.5. 10 pM GST-ybbr-Syt7 C2A 
and ybbr-Syt7 C2A were labelled with the ATTO-565 
via SFP synthase (4/-phosphopantetheinyl transferase). 
Each reaction contained 50 mM tris 7.5, 10 mM MgCl 2 , 
40 mM NaCl, 20 pM ATTO-565 CoA and 1 n M SFP syn¬ 
thase. Reactions were incubated at room temperature for 
30 minutes, then placed at 4°C overnight. Samples were 
dialysed against 1 L of 50 mM HEPES, pH 7.0, 75 mM 
NaCl, 4 mM MgCl 2 and 5% glycerol overnight at 4°C 
then concentrated to 10 fiM. 


Imaging 

All images were acquired using an objective-type to¬ 
tal internal reflection fluorescence microscope (TIRFM). 
The microscope was honre-built around an Olympus 1X71 
body m ina with a 561 nm laser line as excitation 
source. A back-illuminated electron-multiplied charge 
coupled device (EMCCD) camera (Andor iXon DU-888) 
liquid-cooled to -85°C, with an electronic gain of 300 
was used. In order to maintain constant focus during 
the whole imaging time we employed an autofocus sys¬ 
tem (CRISP, Applied Scientific Instrumentation, Eugene, 
OR) in combination with a piezoelectric stage (Z-100, 
Mad City Labs, Madison, WI). Videos were acquired at 
a frame rate of 20 frames/s. 


Image processing and single-particle tracking 

Images were acquired using Andor IQ 2.3 software and 
saved as 16-bit tiff files. Then the images were filtered 
using a Gaussian kernel with a standard deviation of 1.0 
pixel in ImageJ. Single-particle tracking of Atto-C2 and 
Atto-GST-C2 was performed in MATLAB using the U- 
track algorithm developed by Jaqaman et al. [55j under 
thorough manual inspection of detection and tracking. 
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